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Abstract. I review the method of Doppler tomography which translates binary-star 
line profiles taken at a series of orbital phases into a distribution of emission over the 
binary. I begin with a discussion of the basic principles behind Doppler tomography, 
including a comparison of the relative merits of maximum entropy regularisation versus 
filtered back-projection for implementing the inversion. Following this I discuss the 
issue of noise in Doppler images and possible methods for coping with it. Then I 
move on to look at the results of Doppler Tomography applied to cataclysmic variable 
stars. Outstanding successes to date are the discovery of two-arm spiral shocks in 
cataclysmic variable accretion discs and the probing of the stream/magnetospheric 
interaction in magnetic cataclysmic variable stars. Doppler tomography has also told us 
much about the stream/disc interaction in non- magnetic systems and the irradiation of 
the secondary star in all systems. The latter indirectly reveals such effects as shadowing 
by the accretion disc or stream. I discuss all of these and finish with some musings on 
possible future directions for the method. At the end I include a tabulation of Doppler 
maps published in refereed journals. 

1 Introduction 

Many rapidly rotating single and binary stars change little during the course of 
a single rotation or orbit. The spots on single stars can persist for many days, 
while cataclysmic variable stars may stay in outburst for over 100 orbits and in 
quiescence for ten times longer still. However, for the observer, orbital rotation 
can cause considerable variability both in flux and spectra. This arises from 
a combination of changes in aspect angle and visibility, caused by geometrical 
effects, and the rotation of all velocity vectors with the binary orbit. These effects 
are a blessing and a curse: without them we would know considerably less than 
we do about such stars, however the complex variability can be hard to interpret. 

The method of Doppler tomography was developed to unravel the emission 
line variations of cataclysmic variable stars (CVs) [48]. CVs are short period 
binary stars, with orbital periods typically between 1.5 and 10 hours, which are 
beautifully set up to allow us to study accretion. The stellar components of the 
binary, a white dwarf and a low-mass main-sequence star, are faint, and their 
semi-detached configuration means that the geometry is entirely specified by the 
mass ratio and orbital inclination alone. Unfortunately, CVs are far too small 
to be resolved directly - they typically subtend < 10~* seconds of arc at Earth 
- and we can learn nothing of their structure from direct imaging. Instead we 
must turn to more indirect methods. Two key methods are "eclipse mapping". 
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introduced by Horne [39] and reviewed in this volume by Baptista, and "Doppler 
tomography" , the subject of this chapter. Echpse mapping rehes on the geomet- 
rical information contained in eclipse light curves; Doppler Tomography uses the 
velocity information contained in Doppler-shifted light curves. 

In this paper I detail the principles behind Doppler tomography and hope to 
give the reader a full picture of what is now a widely-applied tool in the field. I 
follow the fundamentals with a short section on accounting for stochastic noise, 
which has usually been ignored to date, before finally moving onto a survey of 
results. Although I will attempt to be as self-contained as possible, the subject 
has become too large to cover every application of Doppler tomography and so I 
instead focus upon cataclysmic variable stars. I start this endevour with a potted 
history of the development of Doppler tomography. 

2 History 

Although not restricted in application to CVs, since Doppler tomography was 
developed for CVs and has so far mostly been applied to them (with several hon- 
ourable exceptions covered elsewhere in this volume), for completeness I start by 
describing our standard picture of CVs. Fig. |l| shows a schematic representation 
of our model of non-magnetic CVs with a white dwarf surrounded by a flat disc, 
orbiting a tidally-distorted main-sequence star. In addition a stream of matter 
flows from the main sequence-star and hits the disc in a spot. All of these com- 
ponents, as well as others, have been seen in Doppler images. In eclipsing CVs, 



Fig. 1. A schematic illustration of a cataclysmic variable viewed at four orbital phases. 
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as seen in Fig ^ the main-sequence star occults the accreting regions allowing 
us to locate the chief sources of emission. It was precisely this that led to the 
development of the standard model in the early 1970s [89,72]. In systems with 
compact emission sources, the sharp ingress and egress features during eclipse 
allow the source location to be pin-pointed, and thus, for instance, it is possible 
to determine the binary mass ratio from the observed position of the gas stream 
impact [73,94]. What about more diffuse emission - the disc for instance? An 
obvious possibility is to fit a parameterised model of the disc emission to light- 
curves. However this is plagued by our lack of good a priori models. For instance, 
a symmetric distribution following a power-law in radius is a simple and obvious 
choice for modelling discs, but are discs really symmetric and what if they do 
not follow power-laws? These problems were solved by Keith Home who during 
his thesis work introduced the powerful idea of regularised fitting to the study of 
CVs [39] . Horne modelled the accretion regions with a grid of many independent 
pixels, effectively giving a model of great flexibility. To escape the degeneracy en- 
gendered by such an approach, Horne selected the image of "maximum entropy" 
where the entropy S is given (in the simplest case) by 



where li is the image value assigned to pixel i. This is the method of eclipse 
mapping. 

In eclipse mapping a one-dimensional light curve leads to two-dimensional 
map. Emission lines on the other hand, as I will explain, contain the extra 
dimension of velocity and so can give better constrained images (albeit with some 
disadvantages) . Horne realised that the formation of the line profile from a disc 
was analogous to the formation of medical X-rays used in computed tomography 
to image the human brain. This led to the development of Doppler tomography 
[48], implemented in the first instance in a way analogous to eclipse mapping, 
using maximum entropy regularisation. 

Since that time, and following articles focussing upon the more accessible 
filtered back-projection inversion [40,62], the use of Doppler tomography has 
exploded, with over 100 refereed publications making use of it or containing 
theoretical simulations of Doppler images (Fig. ^). Doppler maps have been 
published for some 16 dwarf novae, 13 AM Her stars, 11 DQ Her stars, 16 nova- 
like variables and 5 black-hole systems, along with other types such as Algols and 
Super-Soft sources. Doppler tomography is approaching the status of a standard 
tool; I now begin discussion of its underlying principles. 




(1) 



Here pi is given by 



(2) 



Pi = 
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Fig. 2. Numbers of refereed publications using Doppler tomograpliy versus year of 
publication. 

3 The principles of Doppler tomography 

Doppler tomography arose from a desire to interpret the emission-line profile 
variations of accretion discs. Past papers on tomography have started by de- 
scribing the formation of the double-peaked profiles from discs [48]. However, 
the existence of a disc is not necessary for Doppler tomography. For instance, 
one of the most spectacular applications of tomography has been to the AM Her 
stars or polars [67] , systems in which the white dwarf has such a strong magnetic 
field that there is no disc. Therefore, in this instance, I will attempt to describe 
tomography from a more general perspective. 

The key to tomography is first to consider a point-like source of emission 
in a binary. Assuming this has a motion parallel to the orbital plane of the 
binary, line emission from such a source will trace out a sinusoid around the 
mean velocity of the system. Assuming one observed such a sinusoid, one could 
associate it with a particular velocity vector in the binary, depending upon its 
phase and amplitude. The trick of tomography is to cope with any number of 
such sinusoids, even when they are so overlapped and blended that one cannot 
distinguish one from another. To understand how this is possible, it is helpful to 
think of profile formation as a projection in the mathematical sense of integrating 
over one dimension of an A^-dimensional space to produce an TV — 1 dimensional 
space. 

3.1 Profile formation by projection 

A given point in the binary can be defined by its spatial position, but, more 
usefully in this case, also by its velocity {Vx,Vy). One must be a little careful to 
define this velocity which is relative to the inertial rather than rotating frame. 
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Inertial frame velocities are always changing as the binary rotates, and so in 
order to define unique values of Vx and Vy, they arc measured at a particular 
orbital phase. Conventionally this is taken to be when the inertial frame lines up 
with the rotating frame. In the case of CVs it is usual to define the a;-axis (in the 
rotating frame) to point from the white dwarf to the mass donor, and the y-axis 
to point in the direction of motion of the mass donor. With this convention, and 
defining orbital phase zero to be the moment when the donor star is closest to 
us, the radial velocity of the point in question at orbital phase ^ is 



where 7 is the mean or systemic velocity of the star. The use of a single value 
of 7 is equivalent to assuming that all motion is parallel to the orbital plane as 
mentioned above. 

With these definitions, an "image" of the system can be defined as the 
strength of emission as a function of velocity, I{Vx,Vy). That is, the flux ob- 
served from the system that comes from the velocity element bounded by to 
Vx + dVx, to + dVy is given by I{Vx,Vy) dVx dVy. I will refer to this as an 
"image in velocity space" or more simply a "velocity-space image" . The relation 
of this image to the conventional image will be discussed below. 

The line flux observed from the system between radial velocities V and V + 
dV at orbital phase (f) can now be obtained by integration over all regions of the 
image that have the correct radial velocity; 



where 5 is a function (of velocity) representing the line profile from any point 
in the image, including instrumental blurring. I assume here that g is the same 
at every point, although it is possible to allow it to vary. The velocity width 
is divided out to obtain a flux density of course, and so the line profile can be 
expressed as 



Ideally g is narrow, best of all a delta function, thus this equation picks out all 
regions of the image close to the line 



This is a straight line inVx,Vy coordinates. Different values of V define a whole 
family of parallel straight lines across the image, with a direction dependent upon 
the orbital phase. With the definition of velocity and orbital phase described 
above, orbital phase corresponds to a collapse in the positive Vy direction, 
phase 0.25 corresponds to the positive Vx direction, phase 0.5 to the negative Vy 
direction etc, with the angle rotating clockwise. Thus the formation of the line 
profile at a particular phase can be thought of as a projection (or collapse) of 



Vr=j — Vx cos 277(1) + Vy sin 27r^, 



(3) 




(4) 




(5) 



V = Vii = J — Vx cos 277(1) + Vy sin 27r^. 
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the image along a direction defined by the orbital phase. Note that if this model 
is correct, two line profiles taken half-an-orbit apart should be mirror images of 
one another. The extent to which this is not the case is one measure of violations 
of the basic assumptions made. 

Phase 0.5 



i 




Veloolly (km/s) 

Fig. 3. A model image and the equivalent profiles formed by projection at angle ap- 
propriate to orbital phases 0.25 (right-most profile) and 0.5 (lower profile). 

Fig. ^ shows a pictorial representation of this process for two projection 
angles. The artificial image has been created with a spot which can be seen to 
project into different parts of the profile at different phases. Tracing back from 
the peaks along the projection directions leads to the position of the original 
spot. This is in essence how line profile information can be used to reconstruct 
an image of the system. 

A series of line profiles at different orbital phases is therefore nothing more 
than a set of projections of the image at different angles. The inversion of pro- 
jections to reconstruct the image is known as "tomography" , the case of medical 
X-ray imaging being perhaps the most famous, although it occurs in many other 
fields too. I now look at the two methods that have been applied in the case of 
Doppler tomography. 

3.2 Inversion Methods 

The mathematics of the inversion of projections dates back to the work of Radon 
in 1917 [60]. If one knows the function (in my notation) /{¥,(/)) for all V and 
(j), a linear transformation ~ the Radon transform ~ can produce the desired 
end product, I{Vx,Vy). In reality, things are not so easy, and we never have 
the luxury of knowing the line profiles at all orbital phases, although one can 
get close in some cases. With the advent of fast computers and the develop- 
ment of medical imaging, interest in the implementation of Radon's transform 
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increased greatly in the 1970s, and one particular method, that of "filtered back- 
projection" found favour [63]. The original paper on Doppler tomography used 
an alternative method inherited from eclipse mapping, that of maximum entropy 
regularisation. In this section I describe these two methods, which are both in 
use today. Each has its pros and cons, which I discuss at the end of the section. 



Filtered back-projection The mathematical inversion of Eq. ^ is detailed in 
the Appendix A. The process can be summarized in the following two steps. 
First the line profiles are filtered in velocity to derive modified profiles, f{V, </>). 
The filter is applied through a Fourier transform, multiplication by |s|/G'(s), 
where G{s) is the Fourier transform over V of g{V) and s is the frequency in 
inverse velocity units, and finally an inverse Fourier transform. It can be applied 
to one spectrum at a time and is a fairly fast process. 
The second step is that of back-projection: 



0.5 



fil — Vx COS 27r0 + Vy sin 27r0, ( 



(6) 



An intuitive understanding of back-projection is extremely useful when trying 
to make sense of Doppler maps. There are two ways of imagining the process 
(see Fig. The first one is perhaps the most obvious from Eq. ^ which implies 
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Fig. 4. Two view of back-projection: on the left are paths of integration through trailed 
spectra (spectra plotted as greyscale images with time running upwards and wavelength 
from left to right). A track close to a sinusoidal component gives a spot in the final 
image. On the right three profiles are smeared back along their original projection 
directions to give a spot. 



that each point in the image can be built by integration along a sinusoidal path 
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through "trailed" spectra (spectra viewed in 2D form with axes of phase versus 
velocity). The particular sinusoid is exactly that which a spot at the particular 
place in the image would produce in the trailed spectrum. This view is illustrated 
on the left of Fig. |. 

However, back-projection is named for another, perhaps more useful, way of 
regarding this operation shown on the right of Fig. ^. In effect, Eq. ^ means that 
the image is built up by smearing each filtered profile along the same direction as 
the original projection which formed it. This way of looking at back-projection 
shows very clearly why small numbers of spectra cause linear artifacts in Doppler 
maps, and one should always be wary of such features. Similarly, any anomalies, 
such as unmasked cosmic rays, dead pixels, flares or unmasked eclipses are liable 
to cause streaks across Doppler maps. 

If the local line profile g{V) is gaussian then so too is G(s), dropping to zero 
at large s. Thus the filter |s|/G'(s) will strongly amplify high frequencies, and the 
image will be corrupted by noise. This may be familiar when it is realised that 
division by G{s) is just the standard (and noise-sensitive) Fourier deconvolution; 
the presence of |s| in this case only exacerbates the problem. One can just remove 
the G{s) term and thus make no attempt to de-convolve the image. Typically 
one goes further still and the filter applied is |s|W(s), where W{s) is a (typically 
gaussian) "window" function to cut off high frequencies and therefore limit the 
propagation of noise into the final image. The penalty for this is that the final 
image is a blurred version of the true image. A similar trade-off will become 
apparent in the maximum entropy inversion which I turn to now. 

Maximum Entropy Inversion In the original paper presenting Doppler to- 
mography, the focus was upon an alternative to filtered back-projection using 
maximum entropy regularisation [48] . This stemmed in part from the earlier de- 
velopment of eclipse mapping but also because I originally developed Doppler 
tomography in spatial coordinates in which the projection becomes one over a 
set of curves rather than straight lines [47]. However, velocity space is nearly 
universally used now, and the linear inversion has become more commonly used 
on the whole. Is another inversion method useful? I think the answer is yes, and 
will discuss why in detail below. First of all, I describe the maximum entropy 
method in some more detail. 

The application of maximum entropy to Doppler tomography is very similar 
to the eclipse mapping case: a grid of pixels spanning velocity space is adjusted to 
achieve a target goodness-of-fit, measured by x^. In general there are an infinite 
number of such images and so the image of maximum entropy is selected. A 
refined form of entropy which measures departures from a "default" image [39] 
is used: 

M 
^=l 
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Here all symbols are as before with the addition of 

where D is the default image. The default image is important in eclipse mapping 
(see Baptista, this volume) but less so in Doppler tomography. Usually a moving 
default is used, computed as a blurred version of the image. This constrains the 
map to be smooth on scales shorter than the blurring, but fixed by the data 
alone on larger scales. For reasonable data, the choice of default appears to have 
little effect. Indeed, it would be my guess that neither does the form of S, and 
that its most important role in this case is to allow a unique solution to be found. 



Relative merits of the two inversions Table |l| compares the Maximum En- 
tropy Method (MEM) and Filtered Back-Projection (FBP) inversion methods; 
plus and minus signs indicate pros and cons respectively. First of all, in the 



Characteristic MEM FBP 

Controlling parameter FWHM of noise filter 

Processor time — (22 sec)^ + (0.7 sec) 

Comparison with data 4-+ 

Consistency of noise level — * + 

Flexibility +4- 

UOO by 100 image, 100 by 50 data, 300 MHz Pentium II. 

Table 1. Comparison of maximum entropy and filtered back-projection for Doppler 
tomography. 



maximum entropy method (MEM) is identified as the equivalent of the window 
filter in the filtered back-projection (FBP). For instance a low forces the 
image to become highly structured in order to fit the data better. Conversely, 
a high allows the image to become smooth and blurred. Next I compare 
processing time, which is perhaps the major disadvantage of MEM. In the ex- 
ample given, MEM took 30 times longer than FBP. Even so, for single images 
the absolute amount of time taken is not large, especially when compared with 
the steps taken to get the data in the first place, although it can become more 
significant when trying to estimate noise, as I will discuss later. Nevertheless, I 
still regard it as a relatively minor disadvantage, and award a single minus. The 
next entry "Comparison with data" refers to the central role that plays in the 
MEM reconstructions which allows one to compare the predicted data directly 
to observations. Filtered back-projection does not try to achieve a good fit to 
the data, leaving one uncertain as to how much better the fit could have been; 
I regard this as a significant disadvantage of the method. 

Propagation of noise into MEM images can be problematical when one is 
comparing, say multiple datasets. The reason is that if one sets the same 
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level for each inversion, in one case this might be easy to reach and a rather 
smooth looking map is the result, while in the next may have a hard time 
reaching the desired level at all, resulting in a noisy map. FBP on the other 
hand seems to give an images of similar appearance for the same window filter. 
It was for this reason that Marsh & Duck [52] used FBP in their tomography 
of the DQ Her star, FO Aqr. It is possible that this problem could be fixed by 
iterating towards a fixed entropy and minimising x^- While this may seem a bit 
ad hoc, it may provide a more accurate representation of how one actually sets 

in practice: it is rather rare to achieve statistically acceptable values of 
for data of reasonably good signal-to-noise ratio, and often is set so that the 
image is neither too smooth nor too corrupted by noise. 

Flexibility, the final entry in Table ^ is another major plus point of MEM. 
For instance, it is very easy to adapt it to the common case of blended lines 
[48,53]. Steeghs (this volume) presents a nice extension of MEM to allow for 
variation of flux with orbital phase. Another more minor example of this is that 
it is easy in MEM to mask out bad data without the need to interpolate. 

For straightforward cases, I think that there is relatively little to choose 
between the two methods, although the speed of the filtered back-projection 
may be advantageous when many maps are being computed. MEM has the edge 
in difficult cases where modifications of the standard model are needed. 

3.3 Noise in Doppler images 

In the penultimate part of this section, I look at the propagation of noise into 
Doppler maps. To an extent Doppler images carry their own uncertainty esti- 
mates in the degree of fluctuation that one sees in the background, and perhaps 
this has motivated the lack of a more rigorous treatment to date. Moreover, it is 
often the I remarked above, that Doppler tomography cannot achieve a 

good fit to the data, and one must assume that systematic errors are dominant. 
However, there is still a need to understand noise, with a common case being the 
question of the reality of a certain feature. Appealing to the level of background 
noise is not always good enough - for example, any map will have a highest 
point, but how is one to judge whether it is significant? 

It has been said that noise cannot be propagated into Doppler maps [84] . It 
is in fact straightforward to do so. However, the important point to understand 
is that noise in Doppler images is correlated. In more detail it is correlated on 
short scales but less so on large scales. This correlation means for example that 
single pixel variances are more-or-less useless in defining the amount of noise 
on a Doppler image. The correlation is positive on short scales and so Doppler 
maps are effectively noisier than an uncorrelated image with the same variance 
per pixel; the difference is significant. 

One of the best ways to appreciate the correlation is to view animations 
running simulated images in movie form. Since this is not possible in print, I 
illustrate the consequence by plotting the scatter in circular apertures in Fig. ^ 
In this figure, although the RMS scatter of the uncorrelated images was adjusted 
to match the MEM images pixel-by-pixel, the scatter at first grows more quickly 



Doppler Tomography 11 




Radius (pixels) 

Fig. 5. The figure shows the RMS scatter in simulated Doppler images measured in 
circular apertures, as a function of the aperture radius for three sets of simulations. 
The solid line shows results of MEM inversions with = 1-0; the dotted line shows 
the FBP results with a filter of FWHM = 1.0 in terms of the Nyquist frequency; the 
dashed line shows what the result would have been if the noise were uncorrelated with 
the same variance per pixel as the MEM simulations. 

with radius in the FBP and MEM images, leading to an RMS about 3 times 
larger at large radii. This illustrates the positive correlation at small separations. 
Ignoring this correlation would lead to a very significant underestimate of the 
true noise. In this case for radii larger than « 4 pixels, all three lines are roughly 
parallel, demonstrating the weak correlation on large scales. The final point is 
the very similar behaviour of the FBP and MEM plots: the two methods lead 
to a very similar propagation of noise into the image. 

The exact pattern of noise is dependent upon the controlling parameter, 
or the filter FWHM. While driving down may force a better fit and higher 
resolution, a penalty is paid in increased noise. Fig. ^ shows this effect. In this 
case as is lowered from 1.1 to 0.9, the noise at the smallest scales increases 
by a factor of two; large scales are almost unaffected. Very similar behaviour is 
seen for different filter widths for filtered back-projection. 

Figures ^ and ^ were made by adding noise to a simulated dataset and then 
reconstructing images. A similar method was employed by Hessman & Hopp 
in their analysis of GD 552 [35]. In practice, the "bootstrap" method [15] is 
preferable. In this technique, real data are used to generate artificial data by 
randomly selecting, with replacement, N new points from N old points. The 
beauty of this method is that it automatically accounts for the true statistics 
of the data and one is not adding extra noise. Consider then an image with a 
feature one wants to characterise. As long as an explicit measurement of the 
feature can be devised, then multiple bootstrap runs can be used to generate an 
uncertainty as well. 
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Radius (pixels) 

Fig. 6. The figure shows the RMS scatter in MEM simulations measured in circular 
apertures for ^ = 0.9 (solid line), 1.0 (dashed), and 1.1 (dotted). 

As an example of how one might use the bootstrap method, consider an image 
showing evidence for a bright-spot, with some question as to the significance of 
the feature (Fig. First one needs a method of measuring the strength of such 
a feature. The method I use here is first to subtract off the symmetric part of the 
image and then to measure the flux in a circular aperture centred on the spot. 
I then use the following procedure to determine the significance of the feature: 

1. Generate a large number 1000) of new datasets from the original data by 
bootstrap resampling (a fast procedure). 

2. Compute maps for each of these in the same manner as for the true map 
(here the CPU penalty of MEM is paid in full). 

3. Subtract the true map from each boot-strapped map to obtain difference 
maps. 

4. Measure the flux of each difference map in the same way as before i.e. by 
subtracting the symmetric part of each image and computing the flux in a 
circular aperture. 

5. Finally, rank the observed flux relative to the fluxes measured from the 
difference maps. 

This procedure generates a set of values showing the stochastic noise level in the 
circular aperture. If carried out for apertures centred on every pixel in the image, 
one can generate a "signiflcance map" which shows the fraction of fluctuation 
values that the observed value exceeds. The result of such a scheme is shown 
in Fig. ^ where an artificial image shown in the top-left led, after the addition 
of gaussian noise, to the trailed spectrum shown in the top-right. The "true 
map" referred to above is displayed in the lower-left and shows evidence for a 
spot in the upper-left quadrant. Carrying out the bootstrap computations with 
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Fig. 7. The panels show an initial test map (upper- left), the equivalent trailed spectrum 
with noise added (upper-right), the reconstructed image (lower- left) and finally the 
lower-right panel shows contours encircling regions of the reconstruction above 99% of 
the 1000 trials (dashed) and 100% of them (solid). 



1000 trials leads to the significance map of the lower-right. This shows that the 
region of the spot is higher than all 1000 of the simulated datasets, whereas no 
other part of the image is. In this case the circular aperture used has a radius of 
200 kms^"'^. and so can be fitted ~ 100 times over into the image. Taking these 
to be independent, then there is ~ 10% chance that one region will exceed all 
1000 trials. More trials would be needed to establish the reality of the feature 
more firmly, but the principle is clear. 

Similar measurements are easily imagined, for instance, one could perhaps 
fit the position of the spot [35], and subsequently obtain uncertainty estimates 
from bootstrapping. All that is necessary is that the measurement is precisely 
defined and applied in the same way to the true and bootstrap maps. 

3.4 Axioms of Doppler Tomography 

Doppler tomography rests on certain approximations to reality that are, at best, 
only partially fulfilled. Violation of these approximations does not mean that the 
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resulting maps are useless, but everyone who carries out Doppler tomography 
should be aware of its limitations. Thus in this final part of this section, I list 
the "axioms" that underly the method: 

1. All points are equally visible at all times. 

2. The flux from any point fixed in the rotating frame is constant. 

3. All motion is parallel to the orbital plane. 

4. All velocity vectors rotate with the binary star. 

5. The intrinsic width of the profile from any point is negligible. 

Exceptions exist to each of these. For instance, emission on the mass donor vio- 
lates the first axiom, while outbursts clearly contravene (2). Doppler tomography 
is an interpretation of the data within a specific model of binary systems and 
only applies inasmuch as the model itself does. 

4 The application of Doppler tomography to CVs 

I now move on to discuss Doppler tomography in the context of cataclysmic 
variable stars. I will mainly focus upon results, but before doing so I consider 
the interpretation of Doppler maps in the case of CVs. 

4.1 Understanding Doppler maps 

On the basis of the standard model presented in section |^ one can easily predict 
the locations of the various components in velocity-space; Fig. ^ shows some 
of the key components of a CV represented in velocity space. The donor star 
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Fig. 8. A schematic of some key components in velocity coordinates. 
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is assumed to co-rotate with the binary, which means that it appears with the 
same shape in velocity as it does in position coordinates, although rotated by 90° 
owing to the relation v = f2 Ar between velocity and position for "solid-body" 
rotation. This reassuring property is somewhat misleading, since the disc, which 
is very definitely not co-rotating with the binary, ends up being turned inside 
out so that the inner disc is at large velocities while the outer disc appears as a 
ring at low velocity. The gas stream is plotted twice: once with its true velocity 
and once with the velocity of the disc along its path; one can also imagine 
intermediate cases. The positions of all these components is fully specified if the 
projected orbital velocities of the two stars, Ki and K2, and the orbital phase 
are known. The overall scale is set by Ki + K2', their ratio, which is the mass 
ratio q = K1/K2 = M2/M1, defines the detailed shape of the stream and Roche 
lobe. The orbital phase sets the orientation of the image, and if it is not known 
the image will be rotated by an unknown amount relative to the "standard" 
orientation shown in Fig. |^. This is not uncommon: for instance if the orbital 
phase is based upon emission line measurements, it is typically delayed by 0.05 
-0.1 cycles with respect to the true ephemeris. This causes an anti-clockwise 
rotation of the image by an equivalent number of turns. The published map of 
LY Hya [78] is a nice example of this phenomenon. 

Although velocity coordinates simplify the picture of line profile formation, 
it is simple enough to invert into position coordinates - indeed this is how I 
originally computed Doppler images [47] . All that is required is a specification of 
the velocity at every point in the system. However, I abandoned this approach 
for two reasons. First, the translation between velocity and position is often not 
known. In fact, perhaps it is never known, given that it is likely that deviations 
from keplerian flow occur. This means that position maps would require recom- 
putation each time system parameters were updated. Second, the same place 
in the system can produce emission at more than one velocity. This is not an 
abstract possibility, but happens in almost every system that has been imaged. 
There are many examples of bright-spot emission from the gas stream while the 
disc at the same location produces emission at a completely different velocity. If 
such data is imaged into position coordinates on the basis of keplerian rotation, a 
spot of emission would be produced at a spurious location in the disc. Sticking to 
velocity coordinates is a reminder of these potential difficulties of interpretation. 
Only in eclipsing systems is there potential for disentangling such effects. 

While we cannot translate the data into position space, there is no difficulty 
in translating any theoretical model into velocity coordinates. Indeed, ideally, 
the theory-data comparison should be made by predicting trailed spectra, doing 
away with the need for Doppler maps altogether. However, Doppler maps still 
have a role in that theoretical models are not good enough to predict all the 
peculiarities of real systems, and comparison is easier in the half-way house of 
velocity space. 

The idea of translating to velocity space also applies to how one should think 
about Doppler maps. Rather than trying to translate features of maps mentally 
from velocity to position coordinates, one should try to think of various com- 
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poncnts and imagine where they would appear in velocity space. The difference 
may seem slight, but it is a significant one. With that said, I now turn to look 
at some results. 



4.2 Doppler imaging results 

There are now a large number of examples of Doppler tomography, covering CVs 
along with other types of binary as well, such as Algols and X-ray transients, the 
latter being very similar to CVs in many ways [51,8]. Rather than spend space 
covering these in detail when the original references do so already, as do other 
contributions in this volume, I have decided to devote this section mainly to 
highlights based upon a literature survey of as many published Doppler images 
as I could find. The results of this survey are tabulated in Tables ^, ^ and || 
contained in Appendix B where I list systems with published maps, which lines 
were mapped, the spectral resolution, the outburst state and some indication of 
the appearance of the maps. 



Spiral shocks The discovery by Steeghs et al. [77] of spiral shocks in the dwarf 
nova IP Peg is perhaps the most significant result from Doppler tomography 
applied to CVs. These shocks appear to be present in all outbursts of IP Peg, 
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Fig. 9. Spiral shocks in IP Peg [27]. 



including (with hindsight!) pre-discovery outburst data [50]; Fig. ^ shows one 
example. There is corroborating evidence from other systems such as SS Cyg 
[76], V347 Pup [82] and EX Dra [42], although none of these are as convincing 
as IP Peg. 

Stream emission in polars Beautiful work by Schwope and others [66,67,30,68,71] 
has revealed the gas stream in the polar class of cataclysmic variables in which 
the white dwarf is magnetically locked to the mass donor star. Although initially 
ballistic, there is evidence for the influence of the field upon the stream and some 
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Fig. 10. Stream-stripping in HU Aqr [30]. The upper panels show the Hell 4686 line 
and Doppler map; the lower panels show simulated data based upon a simple model in 
which gas is pulled of the stream and threaded onto the magnetic field. 



emission from the gas as it hurtles down towards the white dwarf (see Fig. pi]| ). 
Such work has tremendous scope for teaching us about the stream/magnetosphere 
interaction. Changes in the appearance of maps between high and low states have 
been seen in HU Aqr [67] and further observation of differing states should tell 
us how the plasma/field interaction varies with accretion rate. 

Bright-spots Many systems show bright-spots in Doppler maps, classic cases 
being WZ Sge [75] and GP Com [53] . The locations of the spots are interesting. 
In some cases they line up with the stream's velocity [50] , while others are closer 
to the disc's velocity [93]. Still others adopt a position half-way between the 
disc and stream velocities [49,75] (see also Fig ^l]). There is some evidence for 
extended cooling following the spot [75] and there is potential for examining this 
with lines of differing excitation. 

Emission from the secondary star Many Doppler maps show emission from 
the secondary star [50,49,23,57,66]. Some of this may be related to magnetic 
activity of these rapidly rotating stars, but there is little doubt that much of it 
is caused by irradiation by the inner disc. The disc systems are interesting in 
this respect since one can expect the disc to cast a shadow over the equator of 
the mass donor. This seems indeed to be the case [28,57] and perhaps has the 
potential to tell us about the vertical structure of the disc. 
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Fig. 11. The pecuhar bright-spot of U Gem observed in January 2000 (Unda et al in 
prep). 

Missing discs in novalikes Nova-like variable stars are thought to be classic 
examples of steady-state systems. Spectroscopically however they have proved 
hard to understand. A particular peculiarity of these stars is that it is often 
difficult to see any sign of a disc in these systems. Instead the Doppler maps are 
often dominated by a single structure-less blob. This blob tends to be located 
in the lower-left quadrant of the map for the Balmer and Hel lines, but rather 
closer to the expected location of the white dwarf for Hell. Several explanations 
for these phenomena have been proposed but none of them are compelling in my 
view. 

5 The future of Doppler tomography 

The future development of Doppler tomography splits into extensions of the 
method and the acquisition of further datasets. Dealing first of all with the 
latter, it is evident from Tables ^, ^ and ^ that there are several areas where 
improvements are possible. For instance more lines could be covered, especially 
metal lines (e.g. CaH), ultraviolet and infra-red lines. Ha has received relatively 
little attention, but when it has been looked at, often appears peculiar [81,76]. In 
the case of dwarf novae, the discovery of spiral shocks makes extended coverage 
of outbursts of considerable interest. 

Multi-epoch tomography is probably the most gaping hole because it is hard 
to make a coherent picture of the many one-off maps that have been published 
to date. The resolution of maps published to date is poor or moderate in many 
cases, limiting their usefulness. The ultimate limit is set by thermal broadening, 
but we are very far off this in most cases, especially for heavy elements. Pushing 
to high resolution is not trivial because of the concomitant need to shorten the 
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exposure time owing to smearing, but it is thoroughly feasible. Exposures of 
length t can be thought of as blurring the image by a rotation of 



where P is the orbital period. For a feature a speed K from the centre of mass, 
this will match the spectral resolution AV when 



Thus if AV is lowered by a factor of two, the exposure time must also be reduced 
by a factor two. As a concrete example, consider trying to image the mass donor 
in a system where K2 — 400 km s^^. Typical equatorial velocities in CVs are 
« 100 kms""'^, so we may attempt to obtain data with AV = lOkms""'^. We 
would then need A(j) = t/P< 0.004, equivalent to t < 30 s for a system of P = 2 
hours. For many CVs, this will require 8m-class telescopes. 

Various extensions are possible, and have already been developed in some 
cases. Standard Doppler tomography does not treat the geometry of the donor 
star correctly. It is straightforward to fix this [64] . Bobinger et al. [3] describe a 
method for simultaneous Doppler and eclipse mapping of emission lines in which 
a single image is computed to fit both spectra and light curves of the lines, with 
a keplerian velocity field used to translate between position and velocity space. 
It is difficult to evaluate whether the spectra or fiuxes dominate the final maps, 
but it is clear that spectral information does alleviate the degenerate nature 
of eclipse mapping. Of course, the need to assume a particular velocity field 
is a disadvantage. An attempt has been made to avoid this by simultaneously 
adjusting a spatial image and a position-velocity map to fit spectra of eclipsing 
systems [1] . In this method, spectra out of eclipse serve to fix the velocity space 
image as usual, which is then translated to position space through the eclipse 
information. The technique was able to recover a V oc relation from 

spectra of V2051 Oph, but as developed it could not handle the difficult case of 
the same place producing emission at more than one velocity. Finally Steeghs 
(this volume) describes a new modification which allows orbital variability to be 
included in Doppler images. As this is so common, it has considerable potential. 

6 Conclusions 

I have reviewed the principles and practice of the analysis method of Doppler 
tomography which helps the interpretation of the complex line profile variations 
from close binary stars. The key discoveries from the application of this tech- 
nique are the spiral shocks in outbursting dwarf novae and the stream/magnetic 
field interaction in the polar class of cataclysmic variable stars, but Doppler to- 
mography has taught us much about the stream/disk interaction and irradiation 
of the donor star too. 

For the future, efforts need to be made to acquire multi-epoch datasets for 
tomography as these are wholly lacking at present. Following that higher spectral 
and temporal resolution data are needed to exploit tomography to its limit. 
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Appendix A 



In this appendix I show that, as stated in section 3^, Eq. |^ can be inverted 



by application of the filter |s|/G(s) followed by the back-projection of Eq. |6[ I 
define the Fourier transform F{s) of a function /(x), and its inverse by 



/oo 
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— i27TSX 



dx. 



f{x) - / F(s)e'^^^^^ds. 

— oo 

The frequency s here is measured in cycles per unit x. Now take the Fourier 
transform over V of the hne profile equation, |^: 
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Dividing through by G{s), multiplying by |s| and taking the inverse Fourier 
transform gives the filtered line profiles 
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I{V^,Vy) / \s\e-~''^^'^''-^^UsdV.,dVy. (12) 



Finally, back-project these filtered profiles according to Eq. ^, that is compute 
the integral 

^0.5 

~f{VR,4>) dcp, 





where 

Vn — '-f — Vx cos 27r0 -I- Vy sin 2TT(f>. (13) 

Putting dashes on various symbols to avoid confusion later, then the back- 
projection integral becomes 
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The third hne above follows from the second after transforming from polar coor- 
dinates s and (j) to cartesian Sx — s cos 27r(/) and Sy — s sin 27r0, and using Eqs. |3| 
so that 

1 poo />oc poo 



'O <J J —oo J —oo 

(15) 

and then the integrals over Sx and Sy separate to give the two Dirac ^-functions 
of the penultimate line of Eq. |l^ since 

/"OO 

S{x) = / e^''^'^'^^ ds. 



This justifies the assertions of section |3^ 
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Table 2. Doppler maps of CVs and X-ray novae published in refereed journals as of 
September 2000. 
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See end of Tab. W for notes. 



Table 3. Doppler maps of CVs and X-ray novae published in refereed journals as of 
September 2000. 
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H/3, Hel, Hell 


1; 6 


44 




II 




120 


H7, Hel, Hell, 
Call 


1; 3 


54 


BT Mon 


N 


- 


170 


H/3,7, Hell 


2b; 6; 7 


44, 92 


// 


II 


- 


100, blue; 
36, red 


Ha,/3, Hel, 
Hell 


2b; 6 


74 


GQ Mus 


N 


- 


220 


Hell 


1; 3 


16 


CP Pup 


N 




170 


H/3, Hell 


1 


44, 91 


PX And 


NL 


— 


90 


Ha 


1; 2b 


31 






- 


45 


Balmer, Hell 


2b, c; 7 


79 


V1315 Aql 


NL 


- 


65 


H/3, Hell 


2b; 6; 7 


11 








170 


H/3, Hel, Hell 


2b; 6; 7 


44 








115 


H/3, Hell 


2b; 6; 7 


32 


UU Aqr 


NL 




170 


H/3 


1; 2ab 


44, 45 








115 


Ha,/3,7, Hel, 
Hell 


1; 2ab 


38 


V363 Aur 


NL 




170 


H^, Hell 


2bc; 7 


44 


WX Cen 


NL 




130 


H^, Hell 


1; 3; 7 


18 


AC Cnc 


NL 




170 


H/3, Hell 


2bd; 6; 7 


44 


V795 Her 


NL 




60 


H/3 


1 


20 




II 




70 


Ha,/3,7, Hel, 
Hell 


1 


7 


BH Lyn 


NL 


L 


75 


H/3,7,5, Hel 


2b,c; 7 


12 


// 






120 


H^ 


1; 2c 


37 



For codes see Table IJ. 

+These maps were computed on the spin phase of the white dwarf rather than the 
standard orbital phase and thus I have not attempted to describe their features. 
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Table 4. Doppler maps of CVs and X-ray novae published in refereed journals as of 
September 2000. 



Object 


Type 


otate 


Res. 
kms ^ 


Line(s) 


Features 


Ret. 


BP Lyn 


NL 




120 


Ha,/J,7 


1; 2ad 


36 




II 




80 


Ha, Hel 


1; 2bc; 3 


81 


V347 Pup 


NL 




120 


Balmer 


1; 3; 4? 


82 




II 


- 


120 


H/3 


2ab; 6 


19 


LX Scr 


NL 




170 


H/3, Hel, Hell 


2a, b; 6 


44 


SW Sex 


NL 


- 


170 


H/3, Hcl 


1; 2a; 6 


44 




II 


- 


75 


H/3,7,5, Hel, 
Hell 


1; 2a; 6 


14 


VZ Scl 


NL 


- 


170 


H/3 


2ab; 7 


44 


RW Tri 


NL 




170 


H/3, Hel, Hell 


1, 3 


44 




II 


H 


50 


H/3,7 


3 


80 


DWUMa 


NL 


- 


170 


H/3, Hel, Hell 


2b; 6 


44 




II 


L 


75 


Balmer 


3 


13 


UX UMa 


NL 




170 


H/3, Hell 


1; 6 


44 


HU Aqr 


P 


H 


110 


H7, Hell 


3; 5 


66, 30 


AM Her 


P 


H 


130 


NV, SilV, CHI 
(UV) 


2b; 3; 7 


21 


V884 Her 


P 


H 


200 


Ha,/3, Hel, 
Hell 


2ab 


29 


BL Hyi 


P 


H 


160 


Ha, Hel 


2a 


56 


ST LMi 


P 


- 


70 


NaT, Call 


3 


70 






L 


350 


NaT 


3 


41 


V2301 Oph 


P 


H 


80 


Ha,/3,7, Hel, 
Hell 


3; 5? 


71 


VV Pup 


P 


H 


90 


Ha 


3; 5 


17 


MR Ser 


P 




70 


NaL Call 


3 


70 


QQ Vul 


P 


H 


70 


Nal, Mgll, 
Hell, CI 


3 


10 




II 


H 


100 


HcII 


3, 5 


68 


AR UMa 


P 


M 


100 


Ha,/3, Hcl, 
Hell, Mgl, Nal 


3, 5 


65 


RX0719+65 


P 


H 


300 


H/3, Hell 


2ab; 7 


86 


RX1015+09 


P 


H 


? 


Hell 


3; 5 


4 


RX2157+08 


P 


H 


180 


H/3, Hell 


3 


88 



Notes: 

References are amalgamated when they refer to the same data. 



Type codes: BH = black-hole system; DN = dwarf nova; N = old nova; NL = nova-like 
variable; IBWD = interacting binary white dwarf; IP = intermediate polar; P = polar. 

State codes (where relevant): Q = quiescence; O = outburst; SS = stand-still; F = 
flaring; H = high; L = low; M = middle. 

Feature codes: (1) Ring (which may be from a disc), (2) Spot, (3) Secondary star, (4) 
Spiral shocks, (5) Gas stream, (6) Low velocity emission, (7) Little structure or low 
signal-to-noise. 

Type (2) = "spot" does not necessarily imply stream/disc impact, but just refers to the 
appearance of the image. Entries such as 2a refer to the quadrant the spot is located 
in (if the orbital phase is known) . The quadrants start at the upper- left with "a" , and 
then go anti-clockwise from there. A combination such as 2ab implies a spot located 
on the boundary of the upper-left and lower-left quadrants. 

Type (6) = "low- velocity omission" refers to such features as the slingshot prominences 
seen in IP Peg [76] and the emission at low velocities commonly seen in nova-like 
variables. 



